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^ A total of 628 pb~^ of data collected with the ALEPH detector at centre-of-mass energies from 189 
1) to 209 GeV is analysed in the search for gauge mediated SUSY breaking (GMSB) topologies. These 

topologies include two acoplanar photons, non-pointing single photons, acoplanar leptons, large impact 
^ ■ parameter leptons, detached slepton decay vertices, heavy stable charged sleptons and multi-leptons 
• • . plus missing energy final states. No evidence is found for new phenomena, and lower limits on masses 

of supersymmetric particles are derived. A scan of a minimal GMSB parameter space is performed 
^ and lower limits are set for the next-to-lightest supersymmetric particle (NLSP) mass at 54GeV/c'^ 
. and for the mass scale parameter A at lOTeV/c^, independently of the NLSP lifetime. Including the 

results from the neutral Higgs boson searches, a NLSP mass limit of 77 GeV/c^ is obtained and values 

of A up to 16 TeV/c^ are excluded. 
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1 Introduction 



If supersymmetry (SUSY) were an exact symmetry, the new SUSY particles would be 
degenerate in mass with their Standard Model (SM) partners. As no experimental evidence 
has been found to prove the existence of SUSY particles, supersymmetry must be a broken 
symmetry. In Gauge Mediated SUSY Breaking (GMSB) models, supersymmetry is broken 
at a high energy scale in a hidden or "secluded" sector and is then propagated down to the 
visible sector via the SM gauge interactions [0] . The main motivation for GMSB models lies in 
the fact that they can easily cope with the experimental absence of flavour changing neutral 
currents (FCNC). Gauge interactions are flavour blind and the scale at which SUSY breaking 
is mediated is expected to be well below the scale at which flavour symmetry should be broken. 

From the phenomenological point of view the main difference with respect to gravity 
mediated SUSY breaking models (SUGRA) is that in GMSB the lightest supersymmetric 
particle (LSP) is the gravitino (G) which couples very weakly to the other particles. Assuming 
R-parity conservation, SUSY particles are pair produced in e"'"e~ collisions and subsequently 
decay to their SM partner plus gravitinos. Another important characteristic of these models 
is that the next-to-lightest supersymmetric particle (NLSP) is, in general, either the lightest 
neutralino x or the sleptons I. In GMSB models a non-negligible mixing between tl and tr 
states is expected for moderate and large values of tan [3 (the ratio of the vacuum expectation 
values of the two Higgs doublets) or large values of |/i| (the Higgs mixing mass term). If the stau 
mixing is large, the lightest stau fi becomes lighter than the other sleptons, and also possibly 
lighter than the neutralino, being then the only NLSP. 

The lifetime of the NLSP depends on the gravitino mass (or equivalently on the SUSY 
breaking scale \fF which is proportional to it) 0]: 

0.01 /100GeV\^ / \2 
ctnlsp ~ I 777^ ) cm (1) 

where is the bino component of the Xi ^"^^ = 1 for a £ NLSP. When cosmological 
considerations are taken into account, an upper limit is placed on the gravitino mass [Q: 
rriQ < IkeV/c^ (V^ < 2000 TeV). Thus the gravitino mass can range from (9(10~^) eV/c^ 
to IkeV/c^, which practically implies that any NLSP decay length is allowed. For this reason 
topological searches able to identify long-lived or even stable NLSP's have been developed by 
the ALEPH collaboration [g §, §. 

A previous compilation of all GMSB searches carried out by ALEPH exists with data at 
^/s = 189 GeV 0. In this paper the results of the GMSB searches for all data collected at 
\fs up to 209 GeV are summarised. Other LEP and Tevatron experiments have reported their 
results in Refs. i, |, [11. 



The organisation of this paper is the following. A brief description of the ALEPH detector 
is presented in Section |[ In Section ^ the experimental topologies are reviewed and limits on 
sparticle masses are reported. An update on four-lepton final states and the new selection for 
four-lepton final states when sleptons have lifetime are described in Section |3.3| . In Section ^ the 
scan on a minimal set of GMSB parameters is presented. The sensitivity of these parameters 
to the different search exclusions is analysed and lower limits on the NLSP mass and the mass 
scale parameter A are derived. 
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2 The ALEPH detector and data samples 



A detailed description of the ALEPH detector can be found in Ref. |]TT[, and an account of 



its performance as well as a description of the standard analysis algorithms can be found in 
Ref. |]12[. Only a brief overview is given here. 

Charged particle tracks are measured by a silicon vertex detector (VDET), a multiwire 
drift chamber (ITC) and a time projection chamber (TPC). The VDET has a length of 
approximately 40 cm with two concentric layers of silicon wafers at average radii of 6.5 and 
11.3 cm. The ITC consists of eight drift chamber layers of 2 m length between an inner radius 
of 16 cm and an outer radius of 26 cm. The TPC measures up to 21 space points in the radial 
range from 30 cm to 180 cm and has an overall length of 4.4 m. These detectors are immersed 
in an axial magnetic field of 1.5 T and together achieve a transverse momentum resolution 
(t(pt)/pt = 0.0006pT © 0.005 (px in GeV/c). The TPC also provides up to 338 measurements 
of the ionisation energy loss. It is surrounded by the electromagnetic calorimeter (ECAL), which 
covers the angular range |cos^| < 0.98. The ECAL is finely segmented in projective towers 
of approximately 0.9° by 0.9° which are read out in three segments of depth. The energy 
resolution is a{E)/E = 0.18/\/E + 0.009 {E in GeV). The iron return yoke is instrumented 
with streamer tubes acting as a hadron calorimeter (HCAL) and covers polar angles down to 
llOmrad. Surrounding the HCAL are two additional double layers of streamer tubes called 
muon chambers. The luminosity monitors (LCAL and SICAL) extend the calorimetric coverage 
down to polar angles of 34 mrad. 

Using the energy flow algorithm described in Ref. ||12|, the measurements of the tracking 
detectors and the calorimeters are combined into "objects" classified as charged particles, 
photons, and neutral hadrons. A good track is defined as a charged particle track originating 
from the interaction region (with transverse impact parameter |(io| < 1cm and longitudinal 
impact parameter l^^ol < 5 cm), having at least four TPC hits, a transverse momentum greater 
than 200MeV/c and a minimum polar angle of 18.2°. In order to get the correct charged 



multiplicity, photon conversions are reconstructed with a pair- finding algorithm |12|. Electrons 



are identified by comparing the energy deposit in ECAL to the momentum measured in the 
tracking system, by using the shower profile in the electromagnetic calorimeter and by the 
measurement of the specific ionisation in the TPC. The tagging of muons makes use of the hit 
patterns in HCAL and the muon chambers. 

The data samples used in this paper, collected by the ALEPH detector from 1998 to 2000, 
are given in Table |l[ All selections were developed using Monte Carlo techniques. Simulated 
samples corresponding to at least ten times the collected luminosity of all major background 
processes have been generated. A detailed list of the Monte Carlo generators used can be 
found in Refs. [O, IH. Signal samples were simulated with SUSYGEN |1^. The position of the 



Table 1: Average centre-of-mass energy and corresponding luminosities of the analysed data sample. 



Year 


1998 


1999 


2000 


(v/i) (GeV) 


188.6 


191.6 


195.5 


199.5 


201.6 


205.0 


206.7 


/£dt (pb-i) 


173.6 


28.9 


79.9 


87.0 


44.4 


79.5 


134.3 
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most important cuts was determined using the iVgs prescription |TB[ , which corresponds to the 
minimisation of the expected 95% confidence level upper limit on the number of signal events, 
under the hypothesis that no signal is present in the data. 



3 Review of experimental topologies and results 

The nature of the NLSP determines the final state topologies to be studied in GMSB models. 
All the relevant searches according to the NLSP type and lifetime are listed in Table |^. 



Table 2: Final state topologies studied in the different GMSB scenarios. 



NLSP 


Production 


Decay mode 


Decay length 


Expected topology 


X 


e+e" XX 


X-^lG 


^ ^detector 
A ~ ^detector 
^ ^detector 


Acoplanar photons 
Non-pointing photon 
[Indirect search] 


i 


e+e- ^ E 




^detector 
A ~ ^detector 
^ ^detector 


Acoplanar leptons 

Kinks and large impact parameters 

Heavy stable charged particles 


i 


e+e" -> XX 


X^U^ iiG 


^detector 
~ ^detector 
^ ^detector 


Four leptons 

Four leptons with lifetime 
[Not covered here] 



A detailed description of the selections optimised at 189 GeV can be found in Ref . [0] . The 
same selections are applied here with cut values suitably adjusted to take into account the 



beam energy and luminosity increases. Only the four-lepton selections, detailed in Section ^ 
and the acoplanar lepton selection, recently updated in Ref. |[TH], have been modified. 



3.1 Neutralino NLSP 

In the X NLSP scenario, all supersymmetric decay chains will terminate in % — > G7. Searches 
for pair production of neutralinos decaying promptly and neutralinos with intermediate lifetime 
were described in Ref. [0 and updated in Ref. [Q. For short neutralino lifetimes, the resulting 



experimental signature is a pair of energetic acoplanar photons. When the updated analysis 
is applied to the 189-209 GeV data sample, four candidate events are found with 4.9 events 
expected from background processes. For intermediate x lifetimes, one neutralino may decay 
before reaching the electromagnetic calorimeter, while the other decays outside the detector. 
This scenario results in a topology where the only visible energy originates from a single photon 
which does not point to the interaction vertex. Two non-pointing photon events are found in 
the data sample and 1.0 are expected. Systematic errors have been evaluated and included in 
the results as in Ref. [0]. 

For long x NLSP decay lengths (A ^ ^detector) the neutralino becomes invisible and only 
indirect exclusions are possible. The relationship between the x mass and the slepton and 
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chargino masses can be exploited to put indirect limits on the x mass using the ALEPH results 
on slepton [jl3| and chargino fl^ searches performed within the SUGRA framework. 

The combination of all these analyses allows the exclusion of neutralino masses as a function 
of the neutralino decay length. An example for a particular region in the GMSB parameter 
space is shown in Fig. |l|. 




Figure 1: Excluded neutralino mass at 95% confidence level as a function of its lifetime in the 
neutralino NLSP scenario. The excluded areas are obtained by the scan described in Section ^ for 
negative = 1, Mmess = lO^^GeV/c^, tan/3 = 20 and any A. The short and medium lifetime 

cases, when at least one x decays inside the detector, are covered by the acoplanar photons and non- 
pointing photon searches. For long-lived neutralinos the gravity mediated searches for charginos and 
sleptons are used. 



3.2 Slepton NLSP direct decay 

In the case of a slepton NLSP, the pair-production process e"'"e~ U is expected to be the 
main experimental signature. The signal final state topology depends strongly on the slepton 
lifetime. Four different analyses are performed, each corresponding to a specific range of mean 
decay length. These searches were described in Ref. and updated in Ref. 0. 

If the £ has a short decay length, of the order of a few mm or less, the final state topology 
will be a pair of acoplanar leptons and missing energy, carried away by the two gravitinos. 



This final state is studied in gravity mediated models |T3], where a neutralino of almost zero 
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mass plays the role of the gravitino in GMSB. Therefore, the results obtained for a massless 
neutralino are used. 

Sleptons with intermediate lifetime, which decay inside the detector, may show two possible 
signatures: large impact parameter and kinked tracks. If the slepton decays before the TPC, 
between ~1 cm and 40 cm, the slepton track will not be reconstructed and the final lepton 
track will have a large impact parameter. If the slepton decays within the TPC volume, the 
signature is then characterised by a kinked track. Two different selections are therefore applied 
to the intermediate slepton lifetime case. One event is found with 1.1 events expected from SM 
processes. 

Finally, long-lived sleptons can be detected from their anomalous specific ionisation in the 
TPC. The search for heavy stable charged particles selected one candidate event, while 2.3 are 
expected from background processes. 

The effect of systematic uncertainties on kinematic cuts has been studied as in Ref. 
and limits are derived reducing the efficiency by a total systematic error of 5%. When all 
four independent selections are combined, the 95% confidence level lower limits on the right- 
slepton masses, independent of lifetime, are set at 83, 88 and 77GeV/c^ for selectron, smuon 
and stau, respectively. The selectron mass limit is obtained neglecting the t-channel exchange 
contribution to the production cross section. The stau mass limit as a function of lifetime is 
plotted in Fig. ||. 

3.3 Slepton NLSP cascade decays 

In the I NLSP scenario, if the neutralino pair production is kinematically accessible, the process 
e"'"e~ — > XX ^ U ^ d-d-G UG may provide a very distinctive discovery signal. This cascade 
decay can increase the sensitivity to GMSB signatures in some region of the parameter space. 
This process benefits from quite a large cross section (the x is expected to be mainly bino 
and the e^ is expected to be light) and from a clear experimental signature. Four leptons are 
produced in the final state (two could be soft if the x-^ mass difference is small) and in half of 
the cases the two most energetic leptons have the same charge {x are Majorana particles). 

Depending on the flavour of the slepton in x ^ decays, there are six possible topologies, 
labelled ee, flfi, ff, efl, ef and /if in the case of slepton co-NLSP. In the case of stau NLSP 
only the ff topology is relevant. 

3.3.1 Prompt decays 

Searches for the cascade topologies with negligible lifetime have already been performed by 
ALEPH at a centre-of-mass energy of 189 GeV 0. The selection cuts are revisited in this 
paper to improve the signal efficiencies and account for the increase in centre-of-mass energy 
and luminosity. The revised cuts described in Appendix A are also applied to the data collected 
at 189 GeV and improve the signal efficiency up to 10%. 

The main systematic uncertainties on the background and signal expectations come from 
the number of generated events in the simulated samples (up to 4%). A total systematic 
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70/ \ I ■ 

60 / i / \ 

50: / 1 / \ 

10"^^ 10"^^ 10"^° lO"'^ 10"^ 10"^ lO"'^ 

lifetime (s) 

Figure 2: Excluded tr mass at 95% C.L. as a function of its lifetime (shaded area) from direct 
searches. Dashed curves give the limits from the different topologies. The search for acoplanar leptons 
covers the case of small lifetimes, searches for tracks with large impact parameter and for kinks are 
used in the intermediate range, whereas for very large lifetimes a search for heavy stable charged 
particles is performed. The dotted curve gives the expected limit. 

uncertainty of 2% is evaluated for the variables involved in the selection. To derive the results, 
the selection efficiency has been reduced by 4%. 

The efficiency of the cuts on the signal samples was found to be in the range 65-85% for 
mass differences greater than 3 GeV/c^ in the case of final states not involving tau leptons and 
in the range 40-65% for mass differences greater than 5GeV/c^ in the case of final states 
with tau leptons. The selection efficiencies for the six topologies are shown in Fig. 

The numbers of background events expected from Standard Model processes and of events 
observed in the data are given in Table |^. The largest background contributions are from WW 
and Tr(7) events. 

3.3.2 Short and long decays 

A new analysis has been developed for the case of observable slepton decay length, based on the 
experimental signatures of leptons with large impact parameter and track kinks. Searches were 
developed for each of the six channels, consisting of loose cuts on global event variables, which 
cause very little signal rejection, and of more stringent cuts on individual track properties. The 
latter are intended to select the tracks that come from the decay of the long-lived sleptons. For 



6 



Table 3: Expected Standard Model background and selected candidates for the various cascade 
topologies in the case of negligible slepton lifetime. 



Energv 
(GeV) 


ee 


fljl 


e/i 


rr 


er 


fif 


exp 


obs 


exp 


obs 


exp 


obs 


exp 


obs 


exp 


obs 


exp 


obs 


188.6 


1.33 


2 


0.12 


1 


0.98 


1 


5.23 


4 


1.34 


2 


1.77 





191.6 


0.31 


1 


0.02 





0.29 





1.17 


1 


0.33 


1 


0.29 


1 


195.5 


1.05 


2 


0.07 





0.41 





1.84 


7 


0.72 





0.61 





199.5 


0.88 


1 


0.04 





0.69 





2.19 


2 


0.85 


2 


0.86 


1 


201.6 


0.27 





0.03 





0.17 


1 


0.97 


3 


0.30 





0.33 





205.0 


0.51 





0.07 





0.41 





1.96 


1 


0.72 





0.62 


1 


206.7 


0.80 





0.09 





0.70 





3.13 


4 


1.26 


3 


1.24 


2 


Total 


5.15 


6 


0.44 


1 


3.65 


2 


16.49 


22 


5.52 


8 


5.72 


5 



Table 4: Summary of results for the case of observable slepton decay length in cascade decays: 
the numbers of observed candidate events passing at least one topology and of expected background 
events for both decay length selections. 



Energy 


Short decay length 


Long decay length 


(GeV) 


expected 


observed 


expected 


observed 


188.6 


1.39 





0.51 


1 


191.6 


0.24 


2 


0.08 


1 


195.5 


0.66 


1 


0.20 


1 


199.5 


0.73 





0.20 





201.6 


0.38 





0.10 





205.0 


0.69 





0.16 





206.7 


1.16 


2 


0.26 


1 


Total 


5.25 


5 


1.51 


4 



each topology, two independent selections were designed to ensure good sensitivity to the signal 
both for short (~lcm) and long (~lm) decay lengths. Further details are given in Appendix 
B. 

Because the searches in all topologies focus primarily on selecting tracks with large do, there 
is a large overlap in their acceptances, and so candidate events tend to be selected more than 
once. The "efficiency" (the probability for a signal event to be selected at least once) reaches a 
maximum of ~80% for the ff channel and >90% for other channels at a slepton decay length 
of around 10 cm. An efficiency >10% is maintained for slepton decay lengths from ~lmm to 
~10 m for all channels. 

The numbers of background events expected to be selected in at least one topology and the 
corresponding numbers of observed events for each LEP energy and each selection are given in 
Table ^. The efficiency dependence on the slepton decay length is shown in Fig. ^d, for the 
prompt decay and the short and long decay length selections. 
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Figure 3: a) Selection efficiencies for the six different event topologies versus Am = m^-mi in the 
zero slepton lifetime case. The spread of points observed for a given topology is due to different values 
of neutralino and slepton masses, b) Probability for a signal stau-pair event to be selected by at least 
one of the six topological searches versus slepton decay length. The set of curves with higher efficiency 
in the 0.1 cm area corresponds to the prompt decay selection. Those peaking at ~10 cm correspond to 
the short and long decay length selection. Different lines correspond to different points in the (m^,mf ) 
space. 

4 Interpretation of the results in the minimal GMSB 
model 

A scan over the parameters of a minimal GMSB model has been performed to study the impact 
of the different searches. The aim of this scan is to understand which topologies contribute 
to exclude regions in the parameter space and to be able to set a lower limit on the mass of 
the NLSP and on the universal mass scale A, independently of the NLSP lifetime (i.e., for all 
gravitino masses). 

The parameters needed to specify a minimal GMSB model are 

• A, the universal mass scale of SUSY particles; 

• A^5, the number of messenger pairs; 

• Mmess, the common messenger mass scale; 

• tan (3, the ratio between the vacuum expectation values of the two Higgs doublets; 
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Table 5: Minimal set of parameters and their ranges of variation in the scan. 



Parameter 


Lower limit 


Upper limit 




10^ GeV 


10^2 GeV 




10"^ eV 


10^ eV 


A 


103 GeV 




tan/3 


1.5 


40 


N, 


1 


5 


sign(/i) 




+ 



sign(/i), where /i is the higgsino mass parameter; and 

the SUSY breaking scale in the messenger sector, related to the gravitino mass by 



rrir 



2.4 



Planck 



100 TeV 



eV, 



(2) 



where Mpianck = 2.4 x 10^^ GeV/c^ is the reduced Planck mass. 

The ranges of the scan are listed in Table |. The six parameters listed determine the 
properties of supersymmetric particles characteristic of GMSB models. At each point in 
the scan, the ISAJET 7.51 program [|l8l was used to calculate SUSY masses, cross sections, 
branching ratios and lifetimes, then taken into account to evaluate whether a point is excluded 
by any of the searches. In total, over 2.3 million points in the minimal GMSB parameter space 
have been tested. 

In addition to the analyses described in this paper, other searches were used to set exclusion 
areas in the scan: the SUGRA chargino |]l^ and slepton |jl3| searches to cover the case of x 
NLSP with a long-lived neutralino and LEPl results [|^, |20|, used here to exclude very low 
NLSP masses. In addition, for each set of GMSB parameters the Higgs boson masses and 
couplings were computed. The results from Ref. were used to extend the GMSB exclusion 
domain. 



4.1 Lower limit on the NLSP mass 

As reported in Section ^, no significant deviation from the SM expectation was observed. A 
lower limit on the fi mass of 77GeV/c^ is set, independently of its lifetime. This limit is valid 
in the stau NLSP case over the full scan range. It is reduced to 72GeV/c^ in the x NLSP 
scenario. 

The interplay of the different searches in the (m^,mf) plane is shown in Fig. ^. For short 
NLSP lifetimes (Fig. |a) the multi-lepton search is able to exclude up to 92GeV/c^ in 
the slepton NLSP case, extending the acoplanar lepton search. In Fig. ^ the case of long 
NLSP lifetimes is presented. Because they rely on results from indirect constraints, limits in 
the long neutralino lifetime case are less constraining than those obtained with short neutralino 
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ALEPH Vs=189-209GeV ALEPH Vs=189-209GeV 




X mass (GeV/c^) % mass (GeV/c^) 



Figure 4: Regions excluded by the different analyses described in the text at 95% confidence level 
in the (m^,mf) plane, for a) short NLSP lifetimes (m^ < lOeV/c^) and b) long NLSP lifetimes 
(ttIq > IkeV/c^). Points in the dark region are not accessible to the scan. The absolute NLSP mass 
limit is set at 54GeV/c^ in b) by the intersection of chargino and slepton searches. 

lifetime searches. The absolute lower limit on the NLSP mass of 54 GeV/c^, determined by the 
chargino and slepton searches, is visible in Fig. Qo. This point is found at = 1, tan/5 = 3, 
A = 39TeV/c2, M^ess = lO^^GeV/c^ and = lO^eV/c^, where the neutralino is the NLSP 
with the i masses around 96GeV/c^ and all other supersymmetric particles above threshold. 

The impact of the neutral Higgs boson searches on the neutralino and stau mass limits is 
shown in Fig. || as a function of tan/3. The NLSP absolute mass limit is 77GeV/c^ obtained 
for large tan f3 and in the stau NLSP scenario. 

4.2 Lower limit on A 

The parameter A represents the energy scale at which the messenger particles couple to the 
visible sector and hence fixes the universal mass scale of SUSY particles. At the Mmess energy, 
gaugino masses scale like N^A, while scalar masses squared scale like N^A"^. The masses at the 
electroweak scale are calculated by means of the renormalization group equations. Once the 
limit for the NLSP mass has been found, the limit on A as a function of can thus be derived. 

The excluded values for the parameter A as a function of tan/? are shown in Fig. ^ for 
different values of N^,. The absolute limit for A is set at around lOTeV/c^. This limit is set 
at iVs = 5, tan/3 = 1.5, Mmess = lO^^GeV/c^ and a large gravitino mass (stable NLSP). The 
neutralino is the NLSP with a mass of 73GeV/c^, slepton masses are around 76GeV/c^ and 
all other sparticles are above threshold. 
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Figure 5: Lower limits on the masses of a) x ^-nd b) fi as a function of tan /3, for different values of 
A^5, as set by the Higgs boson searches. The shaded area represents the minimum excluded area, for 
any N^, as derived from GMSB searches alone. 



The impact of the neutral Higgs boson searches |21[] is shown in the [A, tan (3) plane in 
Fig. These results strongly constrain the allowed values of A for small tan (3. For example, 
for = 1, A up to 67TeV/c^ and tan/5 up to 6 are excluded. 

The lower limit on A independent of lifetime and tan (3 is shown in Fig. ^a as a function of N^,. 
The absolute limit of A > 10 TeV/c^ can be seen here for = 5. When the Higgs boson search 
results are taken into account, the limit on A increases to 16TeV/c^, for rrit = 175GeV/c^. 
With a top mass of 180GeV/c^, the absolute lower limit on A is set at 15TeV/c^. 

The equation that relates the gravitino mass to the scale of SUSY breaking in the messenger 
sector, niQ = F/ "\/3Mpianck; can be exploited to put an indirect limit on the gravitino mass. The 
universal mass scale must obey A < under the simple assumption of positive messenger 
masses squared |l[]. The lower limit on A can then be converted into a constraint on ^/F 
and therefore provides an indirect limit on the gravitino mass. The dependence of on 
is illustrated in Fig. The lower limit on A of lOTeV/c^ implies a lower limit on rriQ of 
0.024 eV/c^. When the results of Higgs searches are included these limits become 16TeV/c^ 
and 0.061 eV/c^, respectively. 



5 Conclusions 

No evidence for new physics has been found in the search for GMSB topologies in the final 
ALEPH data sample collected at i/i up to 209 GeV. In order to test the impact of the searches 



reported here and in Refs. |13, 17, 21|, a scan over a minimal set of GMSB parameters has been 
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Table 6: NLSP mass limits, as derived from the scan 



NLSP 


mass limit (95% C.L.) 


validity 


X 


92 GeV/c2 
54 GeV/c2 


short X lifetime 
any lifetime 


fi 


77 GeV/c2 


any lifetime 


any 


77 GeV/c2 


Higgs exclusion 



performed. The resulting NLSP mass limits can be read in Table ^. 

The scan also provides a lower limit of 16TeV/c^ on the universal SUSY mass scale A and 
an indirect lower limit on the gravitino mass of 0.061 eV/c^. 
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tanP tanP 

Figure 6: Region excluded at 95% C.L. in the (A, tan/3) plane for a) short, b) fong and c) any NLSP 
hfetime. The impact of the Higgs search is included in d). Values of tan/3 less than 3 are excluded for 
large at any mass parameter A, while tan/3 up to 6 can be excluded for N5 = 1. 
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Figure 7: Exclusions at 95% confidence level (a) for A and (b) for ttIq as a function of N^, derived 
from the minimal GMSB scan (shaded). The unshaded bars represent the excluded region when the 
neutral Higgs boson exclusion is applied with mt = 175 GeV/c^. 
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Appendix A: Four leptons with negligible lifetime 
selection 

The experimental topology consists of two energetic (from the ^ —>■ iG decay) and two soft 
(from the x ^ decay) leptons plus missing energy and momentum. The charge of the two 
most energetic leptons is expected to be the same in 50% of the cases. 

In the analysis the presence of at least two energetic leptons is required, where, in the case 
of tau decays, a lepton can also be a jet with small multiplicity and invariant mass. The same 
muon, electron and tau identification as described in is applied. 

All topologies have a common anti-77 preselection, based on the rejection of events with 
low transverse momentum or with energy deposits at small polar angles, which indicate the 
presence of a scattered electron. The anti^7 cuts include: 

• P±/Vs > 0.075 or ip±/^/s > 0.05 and 10^;,, - 90°| > 15° and - 270°| > 15°); 

• ^dift > 5° or > 15°; 

• EJE,^, < 0.45 and {E^JE,^, < 0.30 or p±^J./^ > 0.03); 

• cosO^i^^ < 0.95 



The cut variables are defined as follows: p± is the transverse momentum of the event, 
and are the azimuthal and polar angles of the missing momentum, 9^^ and ^^catt 
are two angles associated with the 77 kinematic hypothesis described in Ref. [^, i^^h is the 
reconstructed neutral-hadron energy, E^^^ is the total reconstructed energy of the event and 
P±nnh is the transverse momentum of the event evaluated without the neutral hadrons. 



A.l ee, /i/x, e/x selections 

Events with 2, 3 or 4 charged tracks are considered in the case of topologies that do not involve 
tau leptons in the final state. At least three identified electrons (muons) are required in the 
ee (/i/i) selections. In the e/i selection the number of identified leptons must again be at least 
three, but the two most energetic leptons must have different flavours and no more than two 
leptons of the same flavour are allowed. These requirements reject most hadronic background 
decays. After the preselection, the cuts listed in Table ^ are applied to reject planar events and 
improve the 77 suppression. If the two most energetic leptons have different charges, further 
selections are necessary to reduce the remaining WW and ZZ background, and in the case of 
ee, the Bhabha background. 

The cut variables are defined as follows: E^^ is the total energy reconstructed in a 30° cone 
around the beam line, AcopT is the transverse acoplanarity of the two most energetic 



leptons, il and i2. Acol is the angle between il and i2 in space and En is the energy of the 
most energetic lepton in the event. In the ee case, an event is accepted if either the cut on the 
energy of the second most energetic lepton, Ei2, or a veto on isolated photons, 7iso, is satisfied. 
Eit is the energy of the leading track in the event. 
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Table 7: Selection cuts applied to the ee, jljl and e/2 topologies. 



Variable 


ee 


fljl 




AcopT 
Acol 


< 0.12 

< 175° 

< 176° 


< 0.12 

< 175° 


< 0.12 

< 175° 




< 0.37 






if Charge 


L) ^ Charge(£2) 




< 0.39 


< 0.40 


< 0.39 


Enl^s 


< 0.331 

veto i 






7iso 

Acol 




< 172° 








> 0.12 and < 0.37 



A. 2 er, /ir, ff selections 

The selection procedures for the ef and /if topologies are based on the acoplanar leptons 
analysis developed in Ref. where a detailed definition of the tau reconstruction algorithm is 
also given. The selection cuts for mixed topologies are listed in Table |^. 



Table 8: Selection 


cuts applied to the ef 


and /if topologies. 


Variable 


ef 


/if 


A'"(charged) 


= (4,5,6) 


= (4,5,6) 


N{i = /i or e) 


A^(e) > 1 


Ar(/^) > 1 




> 2 


> 2 




> 0.06 and < 0.6 


> 0.06 and < 0.6 


Acol 


< 175° 


< 175° 




< 0.02 


< 0.02 


Thrust 


< 0.96 


< 0.96 




> 0.1 and < 0.8 


> 0.1 and < 0.8 


M(Event - Tl)/^/s 


< 0.54 


< 0.55 


Erl/yfs 


> 0.17 and < 0.38 


> 0.14 and < 0.36 


Er2/V^ 


> 0.01 and < 0.26 


> 0.01 and < 0.30 




> 0.04 and < 0.25 


> 0.02 and < 0.25 




< 0.17 


< 0.18 


En/V^ 


> 0.14 and < 0.36 


> 0.12 and < 0.36 


Ee2/Vs 


< 0.24 


< 0.28 


A^„h(rl) 




< 2 



The definition of the variables used is the following: A^(charged) is the number of 
reconstructed charged tracks, A^(^) is the number of identified leptons {£ = e, /i), M^t is 
the invariant mass of the event, Acol is the acollinearity of the two tau jets, E12 is the total 
energy reconstructed in a 12° cone around the beam line, M^i^g is the missing mass of the event, 
M(Event — rl) is the invariant mass of the event once the most energetic tau jet, rl, has been 
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removed, Eti and are the energies of the two tau jets, Mri and Mr2 are the corresponding 
masses. N^^^lrl) is the number of neutral hadrons reconstructed in the most energetic tau jet. 

The cuts on M^^^ and E-^^ allow further 77 suppression while the cuts on A^(charged), A^(^) 
and A^nh(Tl) reject most of the qq background and hadronic decays. The kinematic cuts on the 
jet variables select events containing tau-like jets. 



For the rr topology the events are clustered into four jets using the Durham algorithm p3 



and each jet is checked for consistency with a tau hypothesis. Events satisfying the cuts listed 
in Table ^ are retained. Two alternative strategies are applied depending on wether the sum 
of the charges of the tracks in the two most-energetic tau jets is zero, Charge(jetl + jet2) = 0, 
or Charge(jetl + jet2) = 1 or 2. 



Table 9: Selection cuts applied to the ff topology. 



Variable 


f f 


-In 1/23 


> 2 and < 8.5 


-lny34 


< 9.5 


-In 7/45 


< 11 


-^^good 


> 2 and < 9 




> 3 


N{r) 


> 3 


N 

^ " cnt 


< 2 




> 0.06 and < 0.5 


^12/yi 


< 0.05 




< 0.09 


AcopT 


< 178° 




< 0.3 




> 0.3 and < 0.9 


Charge(jetl + jet2) 


= 


= 1,2 


Thrust 


< 0.96 


< 0.98 


Acop 


< 172° 


< 171° 




< 0.23 


< 0.26 



The cut variables are defined as follows: 1/23 (1/34, I/45) is the y cut for which the event goes 
from 2 to 3 (3 to 4, 4 to 5) jets. A^good is the number of good tracks, A^^j is the number of 
jets containing at least one charged track, N{t) is the number of identified tau jets, A'cnt is 
the number of charged tracks not associated with any tau jet. -Er^ax is the energy of the most 
energetic tau jet in the event, p^^^^ is the momentum of the most energetic identified lepton (e 
or yu) from the most energetic tau jet (including final state radiation). 

The cuts have been optimised to select events with at least three tau-like jets not lying in 
the same plane. They reject most Standard Model processes leaving an irreducible contribution 
arising from WW and rr(7) events. 
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Appendix B: Four leptons with lifetime selection 



In analogy with the neghgible hfetime analysis, described in Appendix A, six selections 
corresponding to the six different final states have been developed. 

The main requirement for selection is that an event has at least one high-do track. The main 
background to such an object is given by secondary interactions of particles originating from 
SM processes with the detector material (e.g. nuclear interactions, photon conversions, etc). 
The rejection of this background is achieved through detailed analysis of the TPC, ITC and 
VDET hits associated to the reconstructed high-rfo tracks. Three variables are used to select a 
high-(io track: the momentum, the do and Xip (the of the track fit to the interaction point, 
normalised to the number of degrees of freedom). The do cut is not applied if the track \zq\ is 
greater than 8 cm. For each topology two sets of cuts have been applied to cope with the two 
cases of short slepton decay length (rf^ ~ 1cm) and of long slepton decay length {d^ ~ Im). 
These cuts are summarised in Table |TD|. For the ee, jlfi and ejl channels the electron and/or 
muon identification is also applied to the high-rfo track. 

Table 10: The cuts on track momentum (p), do and Xip for each channel under each selection. A 
track must have parameters greater than these values to be considered as a good high-do track in the 
corresponding selection. 



Variable 


Selection 


ee 


/i/i 


TT 


e/i 


er 


fIT 


p (GeV) 


low di 
high d^ 


19.0 
19.8 


20.0 
18.6 


4.7 
6.1 


14.9 
20.0 


10.0 
8.5 


8.1 
4.3 


do (cm) 


low di 
high d^ 


0.17 
0.50 


0.50 
0.50 


0.20 
0.50 


0.20 
0.50 


0.24 
0.49 


0.50 
0.48 


Xip 


low di 
high di 


700 
680 


550 
690 


230 
700 


180 
700 


210 
590 


700 
610 



In addition to the request of having at least one high-(io track, an event is also required to 
pass at least one of six selections based on global event variables. The principal cuts are listed 
in Table 0. 

Most analyses use only tracks with a do less than 2cm and a \zo\ less than 10cm, and 
ignore all others. Some of the variables, indicated as primed, are calculated using only tracks 
with these do and \zo\ conditions, while unprimed variables are calculated using all tracks. The 
definitions of the variables are as follows: Nch is the number of charged tracks, A^tpc is the 
number of charged tracks with at least one TPC hit, i?tot is the total energy of the event, mtot 
is the invariant mass of the event, a is the acoUinearity of the event, C2ch is cosine of the angle 
between the two highest momentum tracks, 5*2^ is a/1 — 0.5(/5^ + /3|) (where jSi and (32 are the 
boosts of the two event hemispheres) and $aco is the acoUinearity of the event. 

The cuts marked with a f are not applied if a parent slepton track has been identified (i.e. a 
mother-daughter relationship has been established between one of the high-c/o tracks and a 
track from the primary interaction point). The cuts on C2ch, S!^/^ and ^'^^^ are only applied in 
the case that there is just one good high-(io track, and that no track or hit has been tagged as 
belonging to its parent slepton. Only one of the cuts grouped with a brace need to be passed. 
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Table 11: Global event variable cuts. 



Variable 


ee, nn, e/i 


ef, flf 


ff 


iVch 


>2 and <20 


>2 






>7 


<9 




Ntpc 




<12 


>2 and <14 


N' 




<10 


<11 






>6 GeV t 


>6 GeV and <0.75^ 






>7.7 GeV t 


>7.7 GeV 


a' 




<178° t 


<178° 


C2ch 
^2/3 


>-0.999 and <0.99 
>0.1 1 
<174°r^ 


>-0.999 and <0.99 
>0.1 -1 
<174°r^ 


>-0.999 and <0.99 
>0.2-i 
<174°r^ 



In addition, cosmic-ray events are suppressed by requiring that the event be within 100 ns of 
the bunch crossing. 
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